Pure spin current from one-photon absorption of linearly polarized light in 

noncentrosymmetric semiconductors 
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We show that one-photon absorption of linearly polarized light should produce pure spin currents 
in noncentrosymmetric semiconductors, including even bulk GaAs. We present 14 x 14 k ■ p model 
calculations of the effect in GaAs, including strain, and pseudopotential calculations of the effect in 
wurtzite CdSe. 

PACS numbers: 72.25.Fe, 72.25.Dc, 72.40.+W 
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There has been a growing interest in the manipulation 
of electron spin in semiconductors 0, Q- Electron 
spin provides a new degree of freedom to be utilized in 
novel devices; indeed, several semiconductor 'spintronics' 
devices have been proposed 0. Some of these require a 
spin-polarized electrical current (SPEC) 0- More exotic 
is a pure spin current (PSC), in which there is no net 
motion of charge; spin up electrons travel in one direction 
while spin down electrons travel in the opposite direction. 
PSCs could result from spin pumpingor Hall effects 00, 
or interference of two o ptic al beams 0, ■ Some of these 
have been observed j^, EH UH • 

In this Letter, we show theoretically that a ballistic 
PSC can be generated in noncentrosymmetric semicon- 
ductors, merely by linear absorption. This feature of 
one-photon absorption of linearly polarized light does not 
seem to have been appreciated previously: It occurs even 
in unstrained bulk GaAs, it can be generated by a single, 
weak continuous wave (CW) laser beam, and it arises 
from the symmetry of the crystal itself. 

We consider a semiconductor with filled valence bands 
and empty conduction bands, uniformly illuminated by 
light with photon energy larger than the direct band gap. 
The electric field inside the semiconductor is E (t) = 
E w exp (— itot) + c.c, where E^ is a slowly varying en- 
velope function. The coupling between the light and 
electrons is treated in the long- wavelength limit, and we 
ignore interactions amongst the electrons. This is a com- 
mon starting point for calculations of interband absorp- 
tion in semiconductors 0] , and we adopt it to show that 
PSCs can be generated even within this relatively simple 
model. We use the semiconductor optical Bloch equa- 
tions (SOBEs) solved perturbatively to first order in the 
field intensity The effect we consider does not rely 
on carrier scattering, but to connect with experiments, 
we include an estimate of the spin separation using a 
simple model of momentum relaxation. 

We indicate below that, within these approximations, 
one-photon absorption of linearly polarized light excites 
a distribution of electrons in k-space that is even in k, no 
matter what the symmetry of the material; no net elec- 
trical current results 0]. Further, electrons excited into 
the conduction band at opposite k points will have oppo- 



site spin polarization, resulting in no net spin injection. 
In noncentrosymmetric crystals however, the spin polar- 
ization injected at a given k need not vanish Thus, 
in such a case, there will be a PSC, since the velocities 
of electrons at opposite k points are opposite. 

By treating the field perturbatively, and assuming fast 
interband dephasing, the SOBE gives for the single par- 
ticle density matrix p nm (t) , 



rnrn 1 / _j 



dt' 



E w {?) 



E* (?) 



(1) 



0(*-f) 



(D nv + D mv ) , 



where pnrn is the initial density matrix, equal to 1 
when ti — m is a valence band and equal to oth- 
erwise, u) nm (k) = uj n (k) — u> m (k), hoj n (k) is the en- 
ergy of the Bloch state |nk), the velocity matrix element 
v cl , (k) = (ck| v vk), 7„ m (k) is the dephasing rate, and 
D nv = 5 (u> — u> nv ) — (i / ir) V (u> — u nv )~ l , where V indi- 
cates the principal part. The dephasing rates approxi- 
mate many-body interactions |l3| . 

The photoexcited electron density is N = 
J2 C k/° cc . k (*)■ Defining N cv (k) as the electron 
density excited from band v to band c at k (i.e. 
N = T,kcv N cvO<-)), we find from Eq. (TJ that 

(d/dt + lcc ) N cv (k) = Mi ni) (k) , where 

iv (k)| 2 S (uj cv (k) - uj) , 



(k) = 



n 2 u 2 v 



IE. 



and V is a normalization volume |l2| . Time reversal 
symmetry of the Bloch states results in Kramers degen- 
eracy, uj n (— k) = LUn (k), where a bar above a band index 
denotes the band with the opposite spin. As well, the 
velocity matrix elements satisfy v cv (— k) = — [v e{j (k)]* 
[r?| . Using these two properties, it follows that when 
the light is linearly polarized, Mi ni) (-k) = (k). 
Consequently, the photocurrent injection rate, given by 
e J2u,c,v I v cc (k) - v vv (k)] ivi" lj) (k), is zero for linearly 
polarized light; hence, any spin current must be a PSC. 

The photoexcited electron spin density is S = 
S c c' k ( c 'k| S |ck) Pcc'.k (t), where S is the spin opera- 
tor. We neglect any spin polarization of the holes, since 
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their spin relaxation times are typically very short (lif . 
Note that for excitation near the band edge hu cc ' is ei- 
ther less than a few meV (if c and c' are equal, or are 
the spin-split lowest conduction bands) or greater than 
a few eV (if c or d is a higher conduction band) [l9j . 
At times t longer than the pulse width t^, the integral 
over t' in Eq. (JIJ will be negligible unless uj cc > < 1/tz,. 
For long pulses, a similar argument can be made but 
with 7 CC ' replacing 1/tz,; however, in this Letter, we fo- 
cus on typical ultrafast experiments 0, 0], for which 
one can neglect spin relaxation and carrier recombina- 
tion occurring on longer timescales. Thus the integral 
over t' allows one to neglect coherences other than those 
between spin-split pairs of bands. We expand the re- 
tained coherences in powers of u> cc i juj and keep only the 
lowest order term. We thus neglect the precession of the 
spins due to the spin-splitting of the bands; this is justi- 
fied since the precession period is long compared to the 
momentum scattering time Q. Writing S = X)k^(k)j 
where S (k) is the electron spin density at k we find 

S l (k) = C jl (k) too E l* (*') E l (*0 < wh ere su- 
perscript indices denote Cartesian components, 



C J '(k) = 



ire 2 1 



£ <c'k|S*|ckK c , (k)v l cv (k) 



v,c,c V / 

x [S (lu cv (k) - u>) + 5 (uj c ' v (k) - cj)] , 



and the prime on the summation indicates a restriction to 
pairs (c, c') for which either c' = c, or c and c' are a quasi- 
degenerate pair. Thus the coherence between bands c 
and c is optically excited and grows with the population, 
as in simpler band models that neglect spin splitting 
If one were to neglect the coherence between spin-split 
bands, the net degree of spin polarization, |S|/ (Nfr/2), 
would only be 10% in GaAs compared to the accepted 
50% with Eq. ©. 

From the time reversal properties of the Bloch states, 
the spin matrix elements satisfy (n, — k S |m, — k) = 
— (n, k| S |m, k)* \vi\. Using this property, Kramers de- 
generacy, and the time reversal property of the velocity 
matrix elements, one finds S (— k) = — S (k) for linear 
polarized excitation. Thus there is no net spin injection 
from linearly polarized light, and if S (k) is nonzero for 
some k, there is a PSC. 

A PSC can be quantified by (If — where Iun is 
the current due to up (down) electrons. A more general 
measure that naturally accounts for a distribution of car- 
rier velocities and spins is to use a spin current density 
pseudotensor K ij = l^ffS 3 ^ 0. 

Phenomenologically, the injection of a spin current 
due to one-photon absorption can be written in terms 
of a material response pseudotensor /i, as K 1 ^ — 
T,i, m I E t (*') E u dt '- The s P in current pseu- 
dotensor /i satisfies the intrinsic symmetry = 
{n l i ml \ , it is further constrained by the symmetry of 



the material, since it must be invariant under the point 
group transformations of the crystal. It vanishes for ma- 
terials with inversion symmetry, but can be nonzero for 
the symmetries appropriate to zincblende and wurtzite 
crystals. 

The microscopic expression for /i is derived from Eq. 
^ , with the same approximations made for Eq. J2J ■ Us- 
ing time reversal properties of the Bloch functions, one 
can then show that /i is real, and can be written as 



k.t? c.c' 



X Rc 



[<c'k| |ck) v l vc , (k) vZ (k) + (I <-> m) 



(3) 



In an ultrafast experiment, a PSC can be measured by 
its displacement of up and down spins. To estimate the 
spin separation, we use the optically injected electron 
distribution as a source term in the Boltzmann trans- 
port equation. By neglecting space-charge effects, which 
are negligible for a PSC, the Boltzmann equation can be 
solved in the relaxation time approximation [ldj . This 
approach neglects scattering during optical excitation, 
but for ultrafast experiments the error should not be too 
large. If one measures the spin with respect to the quan- 
tization direction a, the up and down spin populations 
are separated by an average displacement d (a) . We find 
d l (a) = Y2j (4t/?i) K^aP/N, where r is the momentum 
relaxation time. 

To calculate the PSC for bulk GaAs, we use a k ■ p 
model that diagonalizes the one-electron Hamiltonian 
(including spin-orbit coupling) within a basis set of 14 T 
point states and includes important remote band effects 
[2fJ . The 14 states (counting one for each spin) comprise 
six valence band states (the split-off, heavy, and light 
hole bands), and eight conduction band states (the two 
lowest, which are s-like, and the six next-lowest, which 
are p-like). The model contains 13 parameters chosen to 
fit low-temperature experimental data |20| . 

Strain is included bythe deformation potential method 
of Pikus and Bir 21]. The deformation potentials 
amongst valence and lowest conduction band states are 
well established ■ Between the valence and p-like con- 
duction band states, there are two deformati on p otentials 
for the strain we consider, a cv and b cv [23J, |2J, |25j . We 
use b C v = — 2.3 eV, which is an average of tight binding- 
calculations [23j,|24j, and consistent with experiment [2||. 
The parameter a cv couples to the hydrostatic component 
of strain [24|; neither it nor the deformation potentials 
amongst p-like upper- conduct ion band states affect our 
results. We also neglect the small effect of strain on spin- 
orbit coupling [iil ]. 

It is interesting to compare the GaAs results with a 
material of different symmetry: wurtzite CdSe. Since 
k • p methods are less developed for wurtzite crystals, we 
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FIG. 1: Calculated spin separation in unstrained GaAs (thick 
line), in GaAs under 1% tensile (dotted line) and compres- 
sive (dashed line) strain, and in CdSe (grey line). In GaAs, 
transitions from the split-off band begin to contribute at 
~ 350 meV. The inset shows, for CdSe, the angle 6 between 
d(n) and the c axis; 6 = 45° corresponds to d(n) || E w . 



use a relativistic pseudopotential plane wave code within 
the local density approximation [27j ; the use of a different 
method also illustrates the robustness of the PSC effect. 
Quasiparticle effects are included with a scissor correc- 
tion. The matrix elements of v l S 3 were calculated using 
(c'k| v l & |ck) = J2n v ln 0) (™ k l S j l ck > wi th 16 valence 
and 24 conduction bands included in the sum; nonlocal 
contributions to these matrix elements were neglected. 
The Brillouin zone summation used 2600 irreducible k 
points. 

For GaAs, we present results for light linearly polar- 
ized along [001]. In the standard cubic basis, K xx = 
—K vv = k are the only two nonzero components of the 
spin current. Alternately, using a basis where x || [001], 
i || [110] and y || [110], the only two nonzero compo- 
nents of the spin current are K zy = K yz = k. That is, 
there is a net PSC of electrons with spin component along 
[110] and velocity component along [110], and an equal 
PSC with spin component along [110] and velocity com- 
ponent along [110]. Under strain in the [001] direction, 
the point group symmetry of the crystal is reduced from 
Td to £>2d, but K zy = K yz remain the only nonzero com- 
ponents of the spin current. We have considered 1% biax- 
ial strain; under tensile (compressive) strain, the lattice 
dilates (contracts) by 1% in the [100] and [010] directions 
and contracts (dilates) by 0.93% in the [001] direction, 
with the latter determined from the elastic constants Cn 
and C ia 0. 

For CdSe, we present results for light linearly polarized 
at an angle 7r/4 to the c-axis. Defining n as a vector per- 
pendicular to both the c-axis and the light polarization, 
we restrict our attention to d (n) , which lies in the plane 
perpendicular to n. 



The calculated spin separation distances, |d|, for these 
PSCs are plotted in Fig. 1 as a function of photon energy 
above the band gap, where we have assumed a momen- 
tum relaxation time of 100 fs; note that the band gap 
for each case is different. The inset indicates the direc- 
tion of the PSC in CdSe. For comparison, experiments 
on PSCs generated with harmonically related beams in 
GaAs/AlGaAs quantum wells and ZnSe ^3 measured 
spin separations of 20 nm and 24 nm respectively. 

To illustrate the spin-momentum correlation implied 
by the calculation of d, we plot in Fig. 2 calcu- 
lated angular distributions of the injected carrier den- 
sity N(k) = Yl\k\ N (k) and degree of spin polarization 

s ( k ) = S|k| ^ ( k ) /^( k ) for the excited conduction 
band electrons in GaAs. 

In unstrained GaAs, s maX j the maximum over k of 
|s(k)|/ (h/2), rises from zero at the band edge to a maxi- 
mum of 12% at the photon energy when transitions from 
the split-off valence band become allowed. Under strain, 
Smax is largest near the band edge. In Fig. 2(b), s max is 
15% while in Fig. 2(c) it is 44%. 

It is clear that strain can increase the PSC, especially 
at lower energies where strain-induced splitting of the 
heavy and light hole bands increases the spin polarization 
of the photogenerated electrons. Under tensile (compres- 
sive) strain the light hole band moves to higher (lower) 
energy than the heavy hole band. The calculated split- 
ting between heavy and light hole bands at the T point 
is 84 meV for tensile strain and 69 meV for compressive 
strain. The crystal field splitting between heavy and light 
hole bands in CdSe is 38 meV. 

Note that in Eq. we have assumed that v and S 
commute. This is no longer true when the anomalous ve- 
locity term, %[cr X VV) / (4m 2 c 2 ), is included in v. The 
PSC we have described here does not require that term, 
in contrast to other effects @H^|. In fact, in the k • p 
calculation, the anomalous velocity is neglected. Includ- 
ing /c-dependent spin-orbit coupling between valence and 
lowest conduction bands (with Co — 0.16eV A 29]), and 
the associated anomalous velocity, changes the calculated 
spin separation distances by less than 0.1 nm. 

While the calculations of ballistic PSCs that we have 
presented here are appropriate for bulk semiconductors, 
larger PSCs may be possible in heterostructures, which 
can be prepared with large structural asymmetry. As 
well, we note that since only linear absorption is involved, 
this effect could be studied even with CW beams, and 
in nanostructures, where short spin transport distances 
could still have significant consequences. 

In a material with low enough symmetry, one-photon 
absorption of circularly polarized light can generate a 
ballistic current without an applied voltage |l5j; this so- 
called circular photogalvanic effect (CPGE) is a SPEC 
|30| . It is tempting to understand the effect we present 
here as a superposition of opposite SPECs due to two 
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FIG. 2: (Color online) (a) Calculated k-space distributions of photoexcited electrons in unstrained GaAs with excess photon 
energy 300 meV. The surface is the injected carrier density, N(k), and the vectors affixed to the surface give the degree of spin 
polarization, s(k), of electrons excited with crystal momentum in the direction k. The cross sections through the centers of 
the distributions, along with the projection of spins on the cross sections, have been displaced to aid visualization, (b) and (c) 
are as in (a), but at 20meV excess photon energy and for 1% tensile and compressive strains respectively. The reference spin 
in the lower left corner of each panel is 0.2571/2. 



fields with opposite helicity. However, the CPGE cannot 
occur in crystals with zincblcndc symmetry, although it 
can in wurtzite CdSe j^H and strained GaAs [3^ . It is 
worth noting that in addition to describing the PSC we 
present here, the spin current response pseudotensor in 
Eq. © can also describe the SPEC due to the the CPGE. 

Finally, we point out that a PSC can be gener- 
ated even with unpolarized light. Averaging over po- 
larization directions in a plane, one obtains K lJ — 
J Eldt' J2i, m ( Slm ~ /^' m A where h is the prop- 

agation direction of the light. This PSC is smaller than 
that due to linearly polarized light. 
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